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1 Introduction 



Kowalevski gyrostat is a one parameter integrable extension of the Kowalevski top constructed 
in 1987 > PI] • For the Kowalevski gyrostat the Lax representation with spectral parameter 
was found in framework of a general group-theoretical approach to integrable systems using 
the Lie algebras so(3, 2) or sp(4) j2HEj. For the Kowalevski top the spectral curve generated 
by this Lax matrix differs from the original Kowalevski curve. 

Solutions for the top in terms of the Prym theta-functions were obtained by finite-band 
integration technique in jl], where it was said that the gyrostat can be integrated in a similar 
way. We do not know separated variables and separated equations associated with the Lax 
matrices |lj neither for the top nor for the gyrostat. 

Another algebro-geometric approach to study of the Kowalevski top was developed in jTUl Q]. 
In this method detailed analysis of the level surfaces of constant of motion allows to establish a 
birational isomorphism between the Kowalevski flow and the flows of other integrable systems 
that are linearizable on abelian varieties of the same type. In ^U] the Kowalevski top is related 
with the Neumann system and in [[] with the Schottky-Manakov top on so(4). The known 
Lax matrices and separated variables for these systems give rise to the Lax matrices and the 
separated variables for the Kowalevski top. 

The aim of this paper is to extend the Kowalevski treatment of the top to the gyrostat 
and to construct an isomorphism between the Kowalevski gyrostat and the Clebsch case of the 
motion of a rigid body in ideal fluid. As a byproduct one gets new rational and elliptic Lax 
matrices for the Kowalevski gyrostat together with the corresponding integration procedures. 

2 The Kowalevski gyrostat 

Let two vectors J and x are coordinates on the phase space M. As a Poisson manifold M is 
identified with Euclidean algebra e(3)* with the Lie-Poisson brackets 

\ Ji j Jj } E-ijkJk ; \J% ) "Ej } E-ijk-Ek j {^i i } , (2.1) 

where is the totally skew-symmetric tensor. These brackets have two Casimir functions 

3 3 

A = x 2 = J2^l B = (x-J) = J2^kJk. (2.2) 

fc=l k=l 

Fixing their values one gets a generic symplectic leaf of e(3) 

O ab : {x , J : A = a, B = b} , 
which is a four-dimensional symplectic manifold. 
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H top = -(J! + Jl + 2J*) + cx 1 , ceC. (2.3) 



The Hamilton function for the original Kowalevski top is given by 

1 

This Hamiltonian and additional integral of motion 

K top = f i ■ 6, (2.4) 
are in involution and define a moment map whose fibers are Liouville tori in S a b. Here 

£i = z\ — 2c(x\ + ix 2 ) , & = z\- 2c(x\ - ix 2 ) (2.5) 

and 

z\ = Ji + iJ 2 , z 2 = J\ - iJ 2 . (2.6) 



The Kowalevski gyrostat ^3], (31] is an integrable extension of the corresponding top defined 
by the following constants of motion 

1 

'top 



H = H top -\J 3 = -(J* + Jl + 2Jl-2\J 3 ) + cxi, (2.7) 



K = ^ 2 + A\(^J 3 -X)z 1 z 2 -{zi + z 2 )cx 3 ^ (2.8) 

in involution {H , K} = 0. The gyrostat generalization of the Kowalevski top is essential 
because the corresponding additional terms in the Hamiltonian mimic quantum corrections to 
the top HHJ. 

The equations of motion are given by the customary Euler-Poisson equations 

• dH dH dH , . 

X : J = J X — + xx— , x = xx—, (2.9) 
oJ OX OJ 

where x x z means cross product of two vectors. Equations (|2.9jl can be rewritten as 

x( J x J=P Q dH, P -- 

where 

( J 3 -J 2 \ 

x = 






J 3 


-Ji 


-Js 





Jl 


J2 


-Jl 






J 






-X 


o) 









-x 2 


-x 3 





X\ 


x 2 


—X\ 






We do not show the second commuting flow for brevity. 

Remark 1 [HI]: Additional terms in K are closely related with the third order in momenta 
integral of the Goryachev-Chaplygin gyrostat 

G = 2((J 3 - 2\)ziz 2 + {zi + zi)cx 3 ) 

and expressed in terms of the Poisson brackets of £i, £2 

K = K top + zA{£ 1; 6} - I6\ 2 ziz 2 . (2.10) 



3 



3 The Kowalevski variables 



AcB 4c 2 A- K \ n rr n 7 -2-2 



We will introduce variables s\, 2 for the Kowalevski gyrostat step by step following original 
papers [TH] and [THj where the separation of variables for the top was constructed. 

At first we made a transition from initial variables to new variables £^2 ()2.5|) . (|2.6jl and 
organize four constants of motion in the following matrix identity 

AH AcB \ J Jl cx 3 J 3 \ _ AX f J z - 

CX3J3 c 2 x\ J \ ziz 2 {J z - A) - c(zi + z 2 )x 3 

(ZX+Z2) 2 Z 1 Z 2 (Z 1 + Z 2 ) \ I 6 + 6 £l22+62l \ 

z 1 z 2 (z 1 + z 2 ) z\z\ J \ i x z 2 + ii4 + &4 J 

The second step consists of exclusion two variables £3 and J3 using velocities z% = {H, z{\ 

i z x z 2 + z x z 2 T %{z 2 + zi) , . 

x 3 = , ^3 = 7 - + A. (3.2) 

C Zi-Z 2 [zi-z 2 ) 

Similarity transform U t ( ■ )U of the both sides of (|3.1|) with auxiliary matrix U 

( Zi Z 2 \ ( Z\ -1 

brings us to the following matrix identity for the gyrostat 

+ AiX(zi - z 2 ) 



■ 2 • • 

z l —Z\Z 2 



-Z\Z 2 




(3.4) 



R(z 1 ,z 1 ) i2(zi,z 2 ) 
R(z 1 ,z 2 ) R(z 2 ,z 2 ) 



0. 



Here 



i?(zi, z 2 ) = z\z\ - 2iJ(^ + z \) - AcB( Zl + z 2 ) - Ac 2 A + K. (3.5) 
The diagonal entries of the identity ()3.4)1 allows to express variables 6,2 as 



4iAz fc 4i| - R(z k ,z k ) _ 

Sfc — 7 ^ ) « — 1,Z, 

zi-z 2 [z\-z 2 y 

and one gets integrals H and .ff in terms of biquadratic polynomial R ()3.5|) and two pairs of 
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Lagrangian variables Zi j2 and 

4z 1 Z2 + Rjzx 1 z 2 ) , , 

H = " 2( 2l -^ ' (3 - 6) 

{Az 2 -R{z uZl )){Az 2 -R{z 2 ,z 2 )) 

On the level surface of integrals of motion 

E = {A = a, 5 = 6, # = h, K=k} (3.8) 

relations (|3.6|) and (|3.7|) 

$1,2(21,22,21^2, A B,H,K) | E = 

can be considered as equations of motion determining two dimensional dynamical system. 
Unfortunately variables zi t 2 do not commute {zi, Z2} 7^ 0, so one has to look for more convenient 
parametrization. 

Remark 2 Associated with the fourth degree polynomials R(zk, Zk) (|3.5j) 

R(zk, z k ) = a Zk + Aa\z\ + <oa 2 z\ + Aa^z k + a A , a, e R 
differential equations 

ix 1 i2 (3.9) 



v^(^7^l) \/r(z 2 , z 2 ) ' 

originally appeared in the Euler studies of equation of lemniscate and invariance of the corre- 
sponding elliptic integrals [Hj. In particular Euler proved that equations (|3.9J) have an algebraic 
integral 

£{z!, z 2 , s) = ( Zl - z 2 f s 2 - R(z!, z 2 ) s + W = , (3.10) 
where R(z±, z 2 ) is a mixed biquadratic form similar to ()3.5j) 

R(z x , z 2 ) = a Q z\z\ + 2a 1 z 1 z 2 (z 1 + z 2 ) + 3a 2 (z 2 + zj) + 2a 3 (z x + z 2 ) + o 4 , (3.11) 

and 

w _ R(zi, z 2 f - R(z u z 1 )R{z 2 , z 2 ) 

a( Zi - Z2 y 

In algebro-geometric terms [3U], Euler studied automorphisms (ui,Zi) — > (1*2,22) of the 
algebraic curve of genus one 

C: u 2 = R(z,z), (3.12) 
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which change a sign of the corresponding holomorphic form dz/u — > ±dz/u. Thus every 
algebraic curve of genus one is isomorphic to a complex torus (cubic elliptic curve), which 
equivalent to Jacobian of C. These automorphisms are parameterized by points of a smooth 
elliptic curve 

r: r/ 2 = P 3 (s), P 3 (s)=As 3 + gi s 2 + g 2 s + g 3j (3.13) 

where are functions on initial parameters a , . . . , a 4 . According to Weil |3U], if Ok = 

k = 1,2, are two points of C and = (r)k,Sk), k = 1,2, denote two points of T related by 

1 = N 1 + N 2 and 2 = N x - N 2 then 

dzi dz 2 ds\ dz\ dz 2 ds 2 

— + — = —, — - = — . (3.14) 

m u 2 rji ui u 2 r] 2 

It is infinitesimal version of the Weil interpretation of the Euler results. These results are inde- 
pendent on the choice of affine coordinates (u, z) and (77, s) on the curves C and T, respectively. 

The third step of Kowalevski in JH] is to apply automorphisms of auxiliary elliptic curve 
(|3.12|) given in Remark 2 by introduction her famous variables Si t2 

R(z u z 2 ) ± y/R(z 1 ,z 1 )R(z 2 ,z 2 ) f . 

Sl > 2 = ' ( } 

which are transcendental integrals of the corresponding Euler equations ()3.9j) [6 . 
The variables si 2 are eigenvalues of an auxiliary spectral problem 

R(z!,Zi) R{z 1 ,z 2 ) 
R{z 1 ,z 2 ) R{z 2 ,z 2 ) 

that is naturally extracted from (|3.4|) . Its characteristic polynomial 

£(s) = (zi - z 2 f (s - Sl )(s - s 2 ) . (3.17) 

coincides with the Euler algebraic integral ()3.10|) . For its analysis see, e.g., Golubev 0. 
The matrix of eigenfunctions \I/ of the spectral problem (|3.16J) reads 

/ 1 1 \ 



^ = 2sa 1 (z 1 -z 2 ) 2 q, a 1 =\ , (3.16) 




^/R(zi,zi) \/R(z 2 , z 2 ) 



(3.18) 



\ y/R(zi, zj) ^R(z 2 ,z 2 ) J 



The idea of Kowalevski to pass to new variables S\ t2 is appeared to be very fruitful in her 
treatment of the top. For the Kowalevski gyrostat as well as for the top these variables have 
the following main property: 

Proposition 1 Functions s^ 2 A3. 15}) are Poisson commute {s\, s 2 } = 0. 
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For A = the straightforward proof may be founded in [HJ [Hj. For A ^ this unexpected 
and crucial observation was obtained by direct calculation of the Poisson brackets. It allows to 
suggest that for the gyrostat variables Si$ provide an essential step to separation of variables 
and this stimulated us to write down evolutionary equations and the integrals of motion for it 
in terms of Si )2 and their velocities Si j2 ■ 

Remark 3 In her original paper Kowalevski uses also another preferred coordinate system 
(t, w) in which curve T (j3.13|) has a Weierstrass normal form associated with equation Ps(w) = 
Aw 3 + g 2 w + g 3 . i.e. with gi = (see [T^j, p. 188). These variables Wi = Sj — 2H do not 
Poisson commute. The counterpart of ()3.4j) in variables gives as its off diagonal entries a 
nonphysical identity for w\, w 2 and Wi, w 2 instead of (j3.24jl . 

Function £(zi, z 2 , s) ()3.10|) . (j3.17|) is a quadratic polynomial with respect to any of its three 
arguments z\, z 2 , s. Its partial derivatives with respect to one of variables are discriminants of 
the corresponding quadratic equations. Squares of its partial derivatives with respect to one of 
variables are factorized into functions of the rest two ones 



fd£\ 2 

\da) 



(— =R(z k ,z k )P 3 (s), k = l,2. 



\dz k J 



R(z 1 ,z 1 )R(z 2 ,z 2 ), 

Here polynomial -Ps(s) is given by 

P 3 (s) = As 3 -8Hs 2 + AH 2 s-Ks + Ac 2 A s + Ac 2 B 
Because of complete differential of £(s, Z\, z 2 ) ()3.10|) is zero 

0. 



(3.19) 



(3.20) 



d£ , d£ , d£ , 
— ds + ^—dz x + — dz 2 
OS OZ\ oz 2 



one gets relations between the differentials of the variables of both types 

dsi 2 dzi dz 2 



\/ P 3(Sl,2) y/R(zi,Zx) \ / R{z 2 , Z 2 ) 

which are the Euler equations ()3.14|) . 

In matrix form the relations for velocities look like 



(3.21) 



( h_ \ 

s 2 



~1 



22 



(3.22) 



where we denoted for brevity 



tp h = P 3 (s k ). (3.23) 
Signs at square roots in (j3.21j) . ()3.22j) are compatible with definition of Si,s 2 (|3.15J) and ^ 
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Using (|3.5jl - (|3.22|) we can express integrals of motion H and K ()3.6j) in terms of cubic 
polynomial P%{s), variables si )2 and their velocities s^ 2 

H = ( S 1- S T\- £l±ii, (3.24) 

K ic^tt x / : — — ( Si , S 2 



— = (2# + Sl + s 2 )A 2 -Av/=^n — + — 
4 Vv^i 



- 82) - s lS2 + H 2 . (3.25) 

Here the Hamiltonian H and coefficients of integral K at even powers of gyrostatic parameter 
A are easy calculated using definitions ()3.15j) and ()3.22|) only. For linear in A term in K = 
K 2 \ 2 + KiX + K one gets at first 

Due to inverse of ()3.21|) one converts derivatives d/dz\ j2 to d/dsi i2 



y ^ d . 9 ^ R(z 1 ,z 2 ) 

K x = -4z si- 51 — h s 2 - 



V?i(9si <9s 2 y (^i - z 2 ) 



Minding that from definition ()3.15|) one gets si + s 2 = i?(zi, z 2 )/{zi — z 2 ) 2 and including i into 
square root we obtain finally 



Ki = - i { h \l^r +i2 i^)- (3 - 26) 

In the section |f3 we recover these expressions of H fl3.24j) and K (|3.25jl using relation of the 
Kowalevski gyrostat with the Clebsch system. 

Equations ([3.24)1 . (|3.25j) have the form 

$i, 2 (si, s 2 , h, s 2 , A, B, H, K)\x = 

and depend on the commuting variables si 2, their velocities Si )2 and integrals of motion only. 
Excluding one of the velocities we obtain two equations of fourth order in s& 

(ai - s 2 ) 2 s\ + XV^^h - faifk) 2 + A 2 (si + (2g + " + S2)yfc ) ^ = 0, (3.27) 

where k — 1,2 and is given by 

P k = {2H + Sl + s 2 )X 2 + s 2 + 2Hs + H 2 - — . (3.28) 
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At A = the equations (|3.27|) are reduced to the Kowalevski top equations O EE] 

(-lf(si-S2)4 = V%) , fc = l,2, (3.29) 



which admit integration on S (J3.8|) by Jacobi inversion theorem. Here Ps(s) = P3(s)P 2 (s) is a 
fifth order polynomial, Ps(s) is from ()3.20|) and 

p 2 ( s ) = s 2 + 2Hs + H 2 -^ (3.30) 

is a limiting value of (|3.28J) . P 2 (sfe) = Pk\\=o- 

To construct separation of variables for the gyrostat one needs to substitute Lagrangian 
variables Si, s 2 by momenta 7Ti, 7r 2 conjugated to sx, s 2 and to express integrals of motion H, K 
as functions of Hamiltonian variables. For the Kowalevski top, i.e. at A = 0, momenta 7r 1 ,7r 2 
were extracted from evolution equations (|3.29|) in [S] and Jl] . Integrals of motion as functions 
of s\, s 2 , 7Ti, 7r 2 give rise to separation of variables with separated equations of the form 

Ac 2 b 2 2b 2 
s 2 - AH top Si h H top = 2c 2 (a 2 ) cos (2 v / 2^'vr i ), % = 1, 2 , (3.31) 

where x top = 4H 2 op - K top + 2c 2 a 2 . 

At A ^ transition from velocities Sj, S2 to the corresponding momenta is unknown due to 
complicated form of equations of motion ()3.27j) in si^ variables, thus we cannot claim that si j2 
are separation variables. Nevertheless, below we prove that equations ()3.27|) may be solved in 
quadratures using relation of the Kowalevski gyrostat with the Clebsch system. 



4 The Kowalevski gyrostat and the Clebsch system 

Let two vectors I and p are coordinates on the phase space M.. As a Poisson manifold Ai is 
identified with the algebra e(3)* equipped with brackets 

{h,lj} = e ijk l k , {h,Pj} = e ijk p k , {pi,Pj} = 0. (4.1) 

These brackets respect two Casimir elements 

A=(p,p), B=(p,l). (4.2) 

The following integrable case for the Kirchhoff equations on e(3) was found by Clebsch [S] 

*(p) : l = pxQp, p = pxl. (4.3) 

Here Q is a constant symmetric matrix, det Q ^ . Equations of motion (|4.3J) are generated 
by the brackets (j4.1j) and the Hamilton function 

H = U 2 + 1 -(Qp,p). (4.4) 
9 



The second integral of motion reads as 

/C = (QM)-(Q v p,p), (4.5) 

where Q v stands for adjoint matrix, i.e. cofactor matrix. In our case it reads Q v = (det Q) Q _1 . 
The vector field X (J4.3j) is bi-hamiltonian vector field 

X ( 1 J = V dH = V x dK , (4.6) 



where 



and 



P 




P3 ~P2 

-ps Pl | . (4.7) 
£>2 -pi 

Here I stands for 3x3 unit matrix. Poisson matrices Vo and V\ define two compatible linear 
brackets on M. in a standard bi-hamiltonian formulation. 

Remark 4 If B — the flow ()4.3|) is equivalent to that of the Neumann system with the 
Newton equation of motion 

p = -Qp+((Qp,p)-p 2 )p (4-8) 

describing the motion of a mass point on the sphere p 2 = A under influence of the force — Qp. 

4.1 Mapping of the Kowalevski gyrostat flow onto the Clebsch flow 

The idea of the map Kowalevski top flow ()2.9j) onto the Neumann flow ()4.8|) originally appeared 
in Heine and Horosov (TU1I2S1 to the Kowalevski top and was extended to so(4), so(3, 1) in [To] . 

Let us introduce the following complex vector-functions 

/ Ji Jf + Jf + l J 2 + Jl-1\ „ , A , 



and 



, .cx 3 2cx 3 Ji - J 3 ( J? + J| - 1) .2cx 3 Ji - J 3 (^i 2 + J| + 1)\ , . . , 

Zt °P = ( 2J 2 ' 1 2~J 2 I ' (410) 



such that 



A=(p,p)=a 2 , B=(p,ko P )=0- (4.11) 
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We permuted the first and the second entries in original vectors [10J to make gyrostat formulas 
slightly more symmetric. 

In order to describe mapping of the gyrostat flow (|2.9|) onto the Clebsch flow (|4.3|) we have 
to shift vector ^ op by the rule 

I = Ztop + a^A (p+ik x pY (4.12) 

where k = (1,0,0) is a unit vector. In compare with (j4.11j) scalar product of vectors I and p 
for gyrostat becomes differ from zero 

B=(p,l) = a\. (4.13) 

Adding constraints A = a, B = b to relations (jHHJ), ()4.12|) one gets correspondence M ~ Ai 
of the phase manifolds for the Kowalevski gyrostat and the Clebsch system. Initial Poisson 
structure on M ()2.1j) gives rise to the cubic Poisson brackets { •, -}3 on Ai, for instance 

{Pi,Pk}3 = {Pi(x, J),Pk(x, J)} = e ijk p k (l 2 + i/ 3 )(p2 + ip 3 ) ■ 

We shall not use these induced brackets directly and, therefore, remaining brackets are omitted. 

On the other hand linear Poisson structure on Ai (|4.1|) gives rise to the cubic Poisson 
brackets { ■, -}3 on M. So, M and Ai are multi- Poisson manifolds for which we constructed the 
correspondence M ~ Ai such that their linear brackets map to cubic brackets and vise versa. 

Remark 5 For the Kowalevski top and gyrostat variables p, /top and p, I are coordinates on 
the different spaces A^top and Ai with different brackets (j4.1j) {-, -}top and {■, ■} forming two 
sample of e(3) algebra (|4.1|) . With respect to the top brackets {■, -}top the gyrostat variables 
p, I form the central extension of e(3)top which is contracted to e(3)t p in the limit A — > 0. 

Similar to the top ^U] let us introduce symmetric matrix Q depending on integrals of motion 
of the Kowalevski gyrostat and the Casimir elements on the initial algebra ()2.1|) 



-2 



/ —H —id) icb \ 

— icb — | + c 2 x i (| + c 2 x) 
\ icb i (| + c 2 x) | — c 2 x / 



x = a - K/Ac 2 . (4.14) 



It is easy to prove that this matrix remains constant with respect to dynamics of the Clebsch 
system on Ai generated linear and cubic Poisson structures. 

Proposition 2 Let us identify Ai with M by the map {x, J} — > {p, 1} H4-9^ , H4-l^\ ) such that 
the Casimir elements are equal to 

A = a, B = b, A = a 2 , B = a\, (4.15) 
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If the matrix Q is given by \4-14 ), then 




X = P dH = V H = X. 



(4.17) 



The similar equality holds for the second commuting flows of the Kowalevski gyrostat and the 
Clebsch system. 

The proof is straightforward. 

According to ()4.17|) on the space Ai initial linear Poisson structure Vq and cubic Poisson 
structure induced by Pq generate the same vector field X — X with respect to a common 
integral 7i ~ H (14.17)1 . We can embed ()4.17j) in a standard bi-hamiltonian formulation with 
two functionally different integrals of motion 7i and /C (see (j4.fi j) ). if we extend Poisson space 
Ai — M by additional degree of freedom considering A as an independent dynamical variable. 
Similar extension was used by Sklyanin [2^] when he constructed Lax matrix for the quantum 
Goryachev-Chaplygin gyrostat. 

Thus we arrive at one of the main results of the paper: 

Proposition 3 Solutions of the Clebsch problem give rise to solutions of the Kowalevski gyro- 
stat and vise versa. 

We can get solution of the Kowalevski gyrostat problem using either the Kobb-Kharlamova 
quadratures [201 El HI] , or the Kotter solution [T2| of the Clebsch system in theta functions. 

Recall once more that for A = solutions of the Neumann system was identified with the 
Kowalevski solutions of her problem in |TTH ITH] . 

5 Lax representations 

Equations of motion for the Clebsch system may be expressed in a Lax form 



which automatically exhibits constants of motion as eigenvalues of =Sf and leads to the lineariza- 
tion of the flow on the Jacobi or Prym varieties of the algebraic curve det (Jz?(y) — /il) = 0. 
Here y is an auxiliary variable (spectral parameter). 

There are few different Lax matrices for the Clebsch system associated with two different 
integration procedures. These matrices depend on rational and elliptic matrix-functions of the 
spectral parameter. 



d_ 

If 



■Jf(y) = [Sf(y),*/(y)] 



(5.1) 
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The rational 3x3 Lax matrix for the Clebsch system was found by Perelomov [22] 



J? r (y) = Q + Ly-Ny\ K(y) = V Q + L , (5.2) 

where Q is symmetric matrix (j4.14j) . L is given by ()4.7|) and N = p ® p, N^- = PiPj. The 
corresponding spectral curve is equal to 

n(y,fi) = B 2 y A + ^ 2 + (2H-^trQ)/i-/c)y 2 -det(Q-/iI) = 0. (5.3) 

We do not know separated variables for the Clebsch system associated with this curve. If 
3 = the separated variables are well-known j2I] and may be obtained by various integration 
schemes, for instance, by intersection of two algebraic curves related with Lax matrix J£ r (y) 
|15j . The regular way is provided by the Sklyanin method j2Z] that is appeared to be nontrivial 
in the considered case: 

Proposition 4 In the Neumann case at B = the separated variables U\ % % are poles of the 
corresponding Baker- Axiezer function such that 

££ r (y)^ = /i*, and (a, *) = 1, 

with dynamical normalization a = p = {pi,P2,P3) ■ Canonical variables Uk and their momenta 
p Uk lie on the algebraic curve that gives rise to separated equations. 

The proof consists of direct comparison of the known separated variables [21] with poles 
of the Baker- Axiezer functions. Integration procedure of the generic Clebsch system in these 
variables will be considered in the next section. 

Another Lax matrices for the Clebsch system is related with Kotter's approach [TJj. Let 
Q = diag (ax, a 2 , a^) is a diagonal matrix. Introduce two vectors i(/x) and s(/x) 

t = W(//)Z + W v (//)p, s = W(n)p, W=(^I-Q) 1/2 , (5.4) 

where W(/i) = diag(u>i, w 2 , w 3 ) is diagonal and W v is its adjoint matrix. In a special uni- 
formisation of the spectral parameter /i diagonal entries Wk = — ct^ can be considered as 
basic elliptic functions (see [T7| 12"] l28j). 

The equations of motion 

*0) = s(ji) x t(fi) , (5.5) 
may be rewritten in the Lax form (see |2j) using matrix 

3 3/ v 3 

if e (/i) = ^tfc(/i)o-fe = ^2 f w k l k H 1 2 3 pk J er k , ^e(fi) = y^SfcQtQofc, (5.6) 

k=l k=l ^ k J k=1 

where o~k are the Pauli matrices. 
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The corresponding spectral curve reads 

t 2 (w,h) = w 2 -t 2 (/j,) = w 2 - [A^ 2 + (m-AtxQ)n-K)) - 2iV- det (Q - A*I) = . (5.7) 

The linearization of the flow associated with the curve (j5.7j) and expressions of initial variables 
I and p in theta-functions were done by Kotter using new Lagrangian variables 21,2,-31,2 
which satisfy to nice evolutionary equations and may be considered as candidates for separation 
variables. These variables will be studied in the next section. 

According to |3j let us consider one parametric transformation f ^ : so(4) — > e(3) 

ffj, '■ Pi = Wi(Si - Ti) , U= 12 3 (Si + Ti). (5.8) 

Wi 

where S^Ti are coordinates on so(4) = so(3) © so(3) with the Lie-Poisson brackets 

{Si , Sj } = Eijk Sk , {S^ , Tj } = , {Ti , Tj } = SijkTk . (5.9) 
Remark 6 The inverse transformation f ~ 1 reads as 

*Si = 7; — h + TT—Pi, T i = —^ 1 k~—Pi, (5.10) 

2wiW 2 w 3 ZWi ZW1W2W3 ZWi 

where Wk depends on the parameter /1. The map J" 1 is easy generalized to two parametric 
mapping f~l if we substitute Wk(fi) = v V — a>k and Wk(v) = \fv — in the definition and 
Tfc respectively. 

In Section 7 we show that Kotter used namely this mapping to integrate the Clebsch system. 

The mapping is a twisted Poisson map, which identifies two bi-Hamiltonian manifolds 
e(3) and so(4) such that 

e(3) so(4) e(3) so(4) 

{.,.} {-,.} {.,.} — {.,.} 

instead of for the usual Poisson map. 




{-,,}* {.,.}' {•..}* — {-,-}• 

Here second compatible brackets {. , .}* on e(3) are equal to 

{k,lj}* = e ijk wll kj {h,Pj}* = CijkWjPk, {PhPjY = £ijkh- (5-H) 

The polynomial t 2 (v) is a Casimir function of the corresponding Poisson pencil {. , .}„ = {.,.} — 
u{. , .}*, where v ^ fi. For brevity the second linear brackets {. , .}* on so(4) will be omitted 
because they are completely determined by mapping f ^ ()5.8|) . The similar twisted Poisson map 
related to Steklov integrable cases on e(3) and so(4) is discussed in |29|. 
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The map f ^ identifies the Clebsch system on e(3) with the Schottky-Manakov system on 
so(4) and, according to identify the corresponding 2x2 elliptic Lax matrices. At the 
same time the Schottky-Manakov system on so(4) has another 4x4 Lax matrix with linear 
dependence on spectral parameter. It allows us to construct two-parametric family of Lax 
matrices for the Clebsch system 

_ /z/W v + W v LW- 1 W v p \ _ / Wp 

\ -(W v p) T / \ -(Wp) T z/detW 

(5.12) 

-i In 

where W = [a I — Q) and L is given by ([4.7)1 . This family of Lax matrices leads to a spectral 
surface, rather than a spectral curve. The nature of this surface is discussed in p. 

At \x = the Lax matrix (|5.12|) becomes rational with the following spectral curve 



Ts(y,v) = y 4 - tr Q v uy z + (det Q tr Q v 2 - K) y 2 

- detQ^z/ 2 detQ-2H + trQ^z/y-detQ^detQz/ 2 + i3 2 ) =0 



(5.13) 



The separated variables for the Clebsch system associated with this curve are unknown. 
The corresponding solutions in theta-functions were obtained by algebro-geometric methods in 



At v = the Lax matrix ()5.12|) becomes elliptic matrix, which spectral curve coincides with 
spectral curve ()5.3|) of the rational Lax matrix up to transformation y = By. 



5.1 Algebraic curves associated with the Kowalevski gyrostat 

From relation of the Kowalevski gyrostat with the^Clebsch system established in the section 
flU) one gets three Lax matrices Jf r (y), S£ e {p) and Jzf(/i) for the Kowalevski gyrostat. 

The spectral curve of the rational Lax matrix J£ r {y) looks like 

n(y, n) = X 2 y 4 + (> + (/x + H)X 2 - | ) y 2 - , (5.14) 

where Pz{ix) is given by f)3.30j) . It is a biquadratic function of y with coefficients being quadratic 
and cubic polynomials of \x. At A = it is a famous Kowalevski curve and the corresponding 
variables gi ye rise to the Kowalevski separated variables s\p, (see discussion in the next 
section). 

The 2x2 elliptic matrix ^f e (fi) for the Kowalevski gyrostat has another spectral curve 
r 2 (w,fx) = w 2 - L 2 + + H)\ 2 ~ ^ - X^P^ . 
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Associated with this curve expressions of initial variables x and J in theta-functions may be 
obtained using Kotter formulae and results of the section (J3J). 

At A = the two-dimensional family of 4 x 4 Lax matrix Jf(v, p) (j5.12j) was obtained by 
Adler and van Moerbeke P using directly correspondence of the Kowalevski top flow and the 
Schottky-Manakov flow. For the Kowalevski gyrostat third algebraic curve (|5.13|) is given by 

Ac 2 -±\h- c 2 B 2 ^j (v 2 (^Ac 2 -*pjH- [cBv - \){cBu + A) J (yu - 1) 

Associated with this curve solutions of the gyrostat may be obtained using theta-functions 
expressions for the Clebsch variables from and change of variables from the section 3. 

The fourth algebraic curve associated with gyrostat is due to rational Lax matrices of jl]. 
It is given by 

n(y^) = y 4 -y 2 (2^ - H) - A 2 + ±f) + {» 2 - H) 2 - | + c 2 A - ^ . 



6 Integration of the Clebsch system in elliptic coordi- 
nates 



Minkowski [20] identified the Clebsch system with the Jacobi problem of geodesic motion on 
ellipsoid for which elliptic coordinates ui^ were introduced by Jacobi. In 1895 Kobb started the 
integration procedure in the Euler angles and passed to variables £ = tan(6 l /2), v = tan(0/2), 
which are equivalent to variables u\p [13- 111 1959 Kharlamova [TJ] used directly elliptic 
coordinates U\^. for integration of the second flow of the Clebsch system associated with JC. 

In order to explain the method proposed in ^TJ ^] we reproduce some simple formulae. 
Using equations of motion (|4.3|) and the Casimir elements ()4.2|) we express angular momenta I 
via Lagrangian variables p, p 

I = -j( B P + P x P) ■ 
Then we introduce variables as roots of the following function 

e(/i) = (ji- ui)(// -U2)=ti + I — -j tr Q 1 fi + - . (6.1) 

Substituting tt 1>2 and their velocities Wi )2 into ()4.4|) and ()4.5|) one gets the Hamilton function 
H = T + V 



ui -u 2 ( u x 2 u 2 2 \ B 2 If \ 



(6.2) 
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and the second integral of motion 

JC = ( Ul -u 2 ) = u x * \~ u 2\ ) + — [tTQ-u 1 -u 2 )-u 1 u 2 A. 

\ <Pi ¥2 ) VA V V Vi V P2J A v J 

in terms of variables Ui t2 , their velocities ii\ )2 and the cubic polynomial 

(p k = 4 det (Q - it*!) . 

Below this polynomial will be identified with the cubic polynomial iph ()3.23j) for the Kowalevski 
gyrostat for which we used the same notation. 

Excluding one of the velocities from these equations we obtain two equations of fourth 
degree in each of velocities depending on both variables u\ and u 2 

(A(u 1 -u 2 ) 2 uI + B^-Aip!ip 2 Uk + PkVk) (6-3) 

,2 { a -2 (A 2 (u 1 + u 2 -trQ) + 2AH- B 2 ) Vk \ 2 



V u x -u 2 ) 

Here (3k is a cubic polynomial also depending on u\ and u 2 

(3 k = B 2 ( Ul + u 2 + u k - tr Q) + A (AuI + (2H - Ati Q) Ufc - K) . (6.4) 

Momenta conjugated to iti j2 are introduced by p Ul 2 = J?-, where T is kinetic energy ()6.2|) and 
the Casimir operator B depends on velocities iik- 

Equations (|6.3|) were solved in quadratures in [TTJ H2]. We have to underline only that 
variables {1*1,2, Pm 2 } are no ^ the separated variables in the standard meaning. 

Remark 7 After a suitable rotation 

p = Vp, l = Vl, Q -> Q = V^QV- 1 = diag(ai, a 2 , a 3 ) 

which diagonalize the matrix Q, coordinates u 12 coincide with elliptic coordinates on sphere 
p 2 = A defined by 

( \ t \t \ det (Q - fit) ( p\ p 2 p 2 \ 

e[p) = {ji- ui){n - u 2 ) = + + . 6.5 

A \/jL — ai /1 — a 2 p, — a 3 J 

Remark 8 If B = there are considerably simple equations 

(A{ Ul -u 2 ) 2 u 2 k + 4 Vk p k y = k = 1,2. (6.6) 

In this case both integrals are quadratic polynomials in momenta p Ul2 = ± ^ Ul ~^ Uk and the 
corresponding Neumann system belongs to the Stackel family of integrable systems. Moreover, 
variables {u k ,Pk} are separated variables which lie on the spectral curve ()5.3|) of the rational 
Lax matrix. 
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6.1 The Kowalevski gyrostat in s variables 



Without loss of generality we put a — \p\ — 1 . Inserting p (J4.9j) and I (J4.12j) into the generating 
function e(/i) (j6.1J) of -u-variables one gets 

where S(s) ()3.17|) is the generating functions of the s-variables. Combining this fact with the 
Proposition El we have 

Uk = -Sk — H, u k = -s k . (6.7) 

Here s k = {H,s k }i and ii k = {T~C, u k }2 and {}i,2 means the Poisson brackets (j2.1|) on M and 
the Poisson brackets (j4.1|) on M, respectively. 

Inserting constant of motion (j4.16|) . (|4.15jl and variables (|6.7jl into the equations (|6.2J) we 
arrive at the same expressions of integrals H and K (|H.24jl - (|H.25jl in terms of s-variables and 
polynomials ^ lj2 . However, in contrast with the Clebsch system expression for Hamiltonian H 
()3.24j) contains both integrals H and K via polynomial (p(s). Expressions of H and K in si 2 
and s\ t 2 variables may be obtained by solving equations (jH.27j) . Unfortunately, we cannot use 
very complicated resulting formulae to calculate associated with s variables momenta as in the 
Clebsch case. 



7 The Kotter solution of the Clebsch system 

In 1888 Minkowski proved that the Clebsch flow is isomorphic to geodesic flow on the ellipsoid 
[20J. Then in 1891 Schottky found that the Clebsch flow is isomorphic to integrable motion of 
four-dimensional rigid body, which may be integrated in a special case [25J associated with the 
Neumann system. In 1892 Kotter joined these results together and integrated the Clebsch 
system completely. 

Very brief description of Kotter approach was presented in 0. Below we reproduce the 
essence of Kotter derivation in modern notations. Instead of two vectors t(fi) and s(fi) ()5.4)1 
defining pair of the Lax matrices Jf e (fi>) an d <£^(a0 (|5.6j) we take two vectors t(/i) and t(u) 
defining two samples of the first Lax matrix Jzf e and pass to their linear combinations 

f = a(p, v) tip) + b(ji, u) t(u) = W + (/i, v)l + W+(/i, v)p , 

(7.1) 

V = a (v, v) t(p) ~ Kn> v) t(v) = w -(y" ; v) 1 + W_(/i, v)p , 

depending on two auxiliary functions a, b having zero time derivatives a = b = 0. Here we 
denoted for brevity 

W±(/i, v) = a(/i, u)W(jj.) ± b(fi, u)W(u), W±(/i, v) = a(ji, ^)W(/i) v ± b(fi, z/)W v (z/) . 
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Using inverse relations 

I = -Z- 1 (w_ £ - W+ 77) , V = Z r 1 (w_ e - W + 77) , (7.2) 

where 

Z = 2a(/i, !/)&(//, »/) (w(/i)W» - W(z/)W v (/i)) , 
and the Kirchhoff equations on e(3) ()4.3j) one gets the following evolutionary equations 



fj = (C n 77 + Ci 2 x 77 + (C 21 77 + C 22 £) x £, 
f = (CnV + C12O x f + (C 21 r/ + C 22 x 77, 



(7.3) 



where 



/ C n 


Ci 2 \ 




1 C21 


C22 / 


= ( 



Z- 1 [W(/i) + W(z/)] \ / -W+ w 

Z- 1 [W(/j) - W(u)} J \ -W+ w 

By definition we have CnC 22 — Ci 2 C 2 i = 0. 

Let us postulate the following brackets between £ and rj 

{Vi , Vj } = £ijk(Vk + 7&) , ,€j} = £ijk(£k + 7%) , {& , Ci } = Sijhim + 76) , 

(7.4) 

which coincide with the Lie-Poisson brackets on so(4) ()5.9|) after the following change of vari- 
ables 

q Vi-Ci rp _ Vi + £i 

•Ji — rr~r. r , J- % 



2(1-7)' 1 2(1 + 7)' 
here 7 is arbitrary. 

The equations of motion ()7.3|) are Hamiltonian equations with respect to the brackets (|7.4|) 
and the following quadratic Hamilton function 

^=(C U 77,77) + (Ci 2 £,77)-(C 22 £,0, 

where 

C n C u \ = /1 7 \( C n C 12 

C21 c 22 ) \j 1 ) y C 12 C 22 

So, in fact Kotter constructed two parametric transformation f^ v : so(4) — > e(3) defined by 
relations (|7.1j) that identifies the Clebsch system on e(3) with the Schottky-Manakov system 
on so(4), it means that after mapping f^ u the corresponding Kirchhoff equations on e(3) ()4.3|) 
and on so(4) (|7.3|) coincide. 

If a = 6 then 7 = and the Kotter mapping ()7.1|) is equivalent to the twisted Poisson 
map fn ()5.8|) (see Remark ©). In generic case map f^ u ()7.1|) is the twisted Poisson map 
if a{n,v) and b(u, fi) are numerical functions. In order to prove it we have to check that 
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transformation ()7.1|) gives rise to the second Poisson brackets {.,.}* on e(3) and so(4), 
which are compatible with initial ones. 

According to [23 UZj (see also Remarks (JU) and (JHJ)) equations ()7.3|) may be integrated in 
elliptic coordinates if one of the Casimir elements on so(4) is equal to zero 

(6, V) = , £ 2 + ?? 2 = const. 

However in our case we have 

(£, rj) = a 2 (/i, u)t 2 (fi) - b 2 (fi, u)t\u) , e + r] 2 = 2a 2 (/i, z/)* 2 (/i) + 26 2 (/i, */)*». (7.5) 

Using the two parametric Kotter mapping (j7.1|) the necessary relation (£, rj) = may be 
achieved if one puts \i = £>% and z/ = S&, where s^ are any of the roots of the equation 

t 2 (Sj) = det Sf e (Sj) = , j = 1, 2, 3, 4. 

Of course, these substitutions destroy the Poissonity of the mapping f^y As a sequence, 
substituting n,u = Si,s 2 and /i, z/ = £3,54 conversely into ()7.1|) we can construct two different 
pair of orthogonal vectors 



l/i=S3,V=04 



that satisfy dynamical equations (|7.Hjl and algebraic relations 

^ + ^ 2 = 0, (&*7) = 0, and f + r) 2 = 0, (|, 77) = . (7.6) 
Using relations ()7.2)1 it is easy to prove that vectors £, 77 and 77 are linearly dependent 

SM!::£) (»■ <-> 

where 

U = Z(5 3 ,5 4 )Z- 1 (s 1 ,s 2 )(w_W + - W + W_) _1 (s 3 ,s 4 ); 

Inserting (|7.7jl into (|7.fij) we obtain four algebraic equations, which relate six dynamical variables 
£i and r/j. These equations have non-trivial solutions & = <7i (771, 772, 773) for the special choice of 
the functions a(/i, i/) and 6(/i, z/) (|7.1|) only. 

For instance, Kotter uses the following functions 



where 

V>(/i) = (// - Si) (/j - s 2 ) (/i - s 3 ) (/i - s 4 ) , i//{fi) = 
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In this case relation 

W ± (S 1 ,5 2 )W ± 1 ( S3 ,S 4 ) = W ± (5 1 , S2 )W ± 1 (S 3 ,5 4 ) 

allows us to reduce equations f)7.7|) to more simple form £ = D£ and fj = T) l r), such that four 
algebraic equations (|7.6|) give rise to three independent equations only 

j = l j = l j = l \ 3 ' 

Here d 3 are diagonal elements of the matrix D = W + (s3, s^W+^Si, s 2 ) 

. Wj(s 4 ) + iwj(s 3 ) _ Wj(fx) 



Matrix D obeys the property D*D = —1, where D* = W_(s 3 , s^Wl^Si, s 2 ) is its Hermitian 
conjugate. 

Introduce elliptic coordinates z\^ as the roots of equation 

X(z) z-d{ z-di z-d\ 

where x( z ) = ( z ~ d\){z - d\){z - dj) and A = rjf + r}\ + r/f. From equations ()7.8j) and ()7.9|) 
one gets 



, (fi ~ d 2 j){z2 ~ d 2 j) r]j y/zw ( y/P 5 (*i) y/P&z; 



X'{dfj J d x d 2 d^ zi - z 2 Ui(«i - <i|) ^2(^2-^) 

where -Ps(^) is a fifth order polynomial 

PsO) = z(z - c^d|dl) x{z) = z(z - d?^)(z - d\){z - dj)(z - dl). (7.10) 
By definition ()7.9|) velocities of the coordinates z\^ are given by 

( IVY V 2 x(^) y> TO , -19 



* ^ - d : 

Excluding fjj from equations of motion (|7.3|) one gets 



(-l) J (.ei -* 2 )*j = {dwxzj + du> 2 )<\] P^zj) , j = 1,2, (7.11) 
where G?ra>i 2 are values of the constants of motion, which we reproduce in the Kotter form 



dw x = ^ ^ \/ P ?>( z l) V P 5( Z 2) 



J z\ — z 2 \ z\(z\ — d\)(z\ — d\d\d\) z 2 (z2 — d\)(z2 — d\d\d\) 



(7.12) 



du>2 — 2 ^ ^ 



•izx - z 2 \ {z\ - d\)(zi - djdjdl) (z 2 - d\)(z2 - dfc^d§) 
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(see equation (56) in J7j). Applying the Abel-Jacobi inversion theorem to the equations (j7.11|) 
Kotter then expressed initial variables / and p in theta-functions. 

On the next step we have to rewrite coefficients of Ps{z) ()7.10|) and integrals dw\^ (j7.12J) 
as functions of the initial Clebsch integrals and have to prove that the Kotter variables zip 
(j7.9J) are in involution with respect to initial Poisson brackets for the Clebsch system. It 
is not easy task since {f]i,r]j} ^ and {dk,r]j} ^ 0. We suppose that or commute or 
their commutativity may be restored by using another functions a(/i, v) and b(fi, v) in (j7.1|) 
because the Kotter solution in theta-functions was reproduced in framework of the finite-band 
integration technique [2]. Recall that in [HI] Kowalevski used non-commutative variables too 
(see Remark (JHJ). 

On the other hand, the matrix Jzf e (/i) (|5.fi|) belongs to the family of the Lax matrices for 
elliptic or XYZ Gaudin magnet. The separated variables for the generic Gaudin model were 
constructed in [2E] in classical and quantum mechanics. So, we can construct the separated 
variables for the Clebsch system and, therefore, for the Kowalevski gyrostat using the Sklyanin 
method, which is free from the difficulties of the Kotter approach. 

8 Conclusion 

Our treatment shows that including gyrostatic term into the Hamiltonian of the Kowalevski top 
we arrive at the essentially more complicate dynamic system. Recall once more that quantum 
corrections to the Kowalevski top looks similar to the gyrostatic term. 

The main result of this paper is that we establish one to one correspondence between the 
Kowalevski gyrostat and the Clebsch system. It allows us to construct solution of the gyrostat 
problem using various known solutions of the Clebsch model. The separation of variables for 
the Clebsch system that is unknown now becomes a matter of a primary importance. 

The similar correspondence may be obtained for the Kowalevski gyrostat on so(4) algebra, 
which is equivalent to the generalized Kowalevski gyrostat on e(3) ^Hj. For the so(4) top initial 
vector Z so ( 4 ) was obtained in [T^j. In the gyrostat case we have to substitute it by the rule (j4.12j) . 

Authors thank V V Sokolov for revival our interest to the Kowalevski gyrostat and for useful 
discussions. The research of AVT was partially supported by RFBR grant 02-01-00888. 
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